On the clustering of finite-size particles in turbulence 
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We investigate experimentally the spatial distributions of heavy and neutrally buoyant particles 
of finite size in a fully turbulent flow. As their Stokes number [i.e. ratio of the particle viscous 
relaxation time to a typical flow time scale) is close to unity, one may expect both classes of 
particles to aggregate in specific flow regions. This is not observed. Using a Voronoi analysis we 
show that neutrally buoyant particles sample turbulence homogeneously, whereas heavy particles 
do cluster. These results show that several dimensionless numbers are needed in the modeling (and 
understanding) of the behavior of particles entrained by turbulent motions. 
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Turbulent flows laden with particles are widely 
found in industry and nature; their study is there- 
fore of great interest and holds many fundamental 
aspects, issues and limits still to be explored. One 
striking feature of these flows is the trend for the 
particles to concentrate in specific regions of the 
carrier flow. This has been observed and inves- 
tigated for a long time both in experiments and 
simulations, and it is still widely studied (see the 
review paper [1] and references therein). The fo- 
cus is usually put on small (namely much smaller 
than the dissipation scale -q of the flow) and heavy 
particles (with a high density compared to that 
of the fluid), especially in numerical studies. Be- 
cause of their high speciflc density, the dynamics 
of such small and heavy inertial particles deviates 
from that of the carrier flow. Clustering phenom- 
ena are then one of the many manifestations of 
these inertial effects, generally attributed to the 
centrifugal expulsion of heavy particles from turbu- 
lent vortices, and more recently to a sticking effect 
of zero-acceleration points of the carrier flow [51. 
Other studies indicate that light particles exhibit 
the same trend to cluster but with different clus- 
ter geometries [3j H] . Finally tracers (ought to be 
neutrally buoyant and much smaller than r]) are 
used as seeds to characterize the flow dynamics 
and should not cluster. As small and heavy par- 
ticles, finite-size heavy particles have been found 
to cluster [F. However, the case of finite-size neu- 
trally buoyant particles (with a diameter signifi- 
cantly larger than 77) has never been treated to our 
knowledge in the context of preferential concentra- 
tion. Such particles are known experimentally \<o- 
[S] and numerically 0110], to differ from tracers, 
with in particular different acceleration statistics. 
But existing studies have focused on the dynamics 
of isolated particles, not on the spatial structura- 
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tion of laden flows. Whether they cluster or not 
remains an open question. 

Particles interacting with a turbulent flow are 
commonly characterized by their Stokes number, 
that is the ratio between the particle viscous re- 
laxation time Tp and a typical time scale of the 
flow. Dealing with flnite-size particles, we use the 
same deflnition as in Ref . |1H IT^ , using the eddy- 
turnover time at the scale of the particle, tj_, as 
the time scale of the flow, and a corrective factor 
ftp based on the Reynolds number at particle scale: 
St = Tp/Td = /004/3(l + 2r)/36, where T = pp/pf 
is the particle to fluid density ratio and (j) = d/rj is 
the particle diameter normalized by the dissipation 
scale (note that our conclusions remain the same 
using a point-particle definition of the Stokes num- 
ber). This dimensionless number is often used as 
the key parameter to characterize particle dynam- 
ics in turbulence, using simple Stokesian models 
where the dominant force acting on the particle is 
taken as the drag due to the difference between the 
particle velocity and that of the fluid. These mod- 
els predict preferential concentration of particles 
with non- vanishing Stokes number, with a maxi- 
mal segregation for St around unity [Ol [T3] • This 
behavior is confirmed, at least qualitatively, in ex- 
periments with small and heavy particles |15] . In 
the present study we address the case of finite-size 
particles (both neutrally buoyant and heavier than 
the fiuid) and investigate the particle concentra- 
tion field function of their Stokes number, in 
a homogeneous and isotropic turbulent flow. First, 
we describe the experimental setup and the data 
processing used. Then, we describe the results on 
spatial structuration for finite-size particles. We 
finish with a brief discussion and conclusions. 

In order to study behaviors of both neutrally 
buoyant and heavy particles, we use water as the 
carrier fiuid. The turbulence is generated in the 
Lagrangian Exploration Module (LEM, see Fig.[l]), 
whose characteristics are described in detail in 
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Table I: Turbulence characteristics, /imp: rotation frequency of the 12 impellers; u': fluctuation velocity of the 
flow; e: energy dissipation rate; r) = [y^ je)^^'^ and = {v/eY^'^: Kolmogorov length and time scales of the flow; 
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Figure 1: (color onhne) (a) CAD drawing of the 
LEM (b) Schematic upper view of the setup. 



Ref. US] and summed up in Table |ll The LEM 
produces turbulence in a closed water flow forced 
by 12 impellers evenly distributed on the faces of 
an icosahedral vessel. Here, all impellers rotate si- 
multaneously at the same constant frequency /imp 
which can be varied from 2 Hz up to 12 Hz; op- 
posing impellers form counter-rotating pairs. This 
produces a very homogeneous and isotropic turbu- 
lence in a significant central region of the device 
of the order of (15 cm)'^, with the integral length 
of the order of 5 cm [TS] . It also permits to obtain 
high Reynolds number turbulent flow with mean 
velocities much weaker than the fluctuations near 
the center of the apparatus. Acquisitions are per- 
formed using 8 bit digital imaging at a resolution 
corresponding to a visualization window of the or- 
der of 16 X 12 cm in the center of the LEM. The vi- 
sualization zone is illuminated by a 100 W Nd:YAG 
pulsed laser synchronized with the camera, creat- 
ing a green light sheet with millimetric thickness. 
The camera is equipped with a macro lens. A 
Scheimpflug mount compensates for the depth of 
field effects resulting from the angle between the 
camera and the laser sheet (see Fig. [ij. Images 
are recorded at a low repetition rate of 2.5 Hz suf- 



ficient to address particle spatial distribution as we 
ensemble average over 2000 independent fiow real- 
izations for each experiment (particle dynamics is 
not addressed here). 

We explore the behavior of two kinds of finite- 
size particles, neutrally buoyant and heavy ones. 
Regarding the neutrally buoyant particles, we use 
expanded polystyrene particles whose density has 
been irreversibly adjusted (prior to actual exper- 
iments) by moderate heating so that the density 
ratio to water is 1 < Fn < 1.015. Depending on 
their size, the particles can be regarded as tracers 
following the carrier fiow dynamics, or as finite-size 
particles whose dynamics departs from that of the 
fiow. The limit for tracer behavior is known to be 
c?n ~ 5?7, whereas finite-size effects appear for cZn > 
57y [3 m [13 • In our case the particle diameter is 
c?n = 700 |xm ± 20 |xm and ranges from 4.5?? to 17r; 
(from the lowest to the highest investigated R\). 
Hence, our particles transit from tracers to finite- 
size particles as /imp increases. As for heavy parti- 
cles, we use slightly poly-dispersed sieved glass par- 
ticles with diameters 225 |xm ±25 |J.m and a density 
ratio Fh = 2.5, making them inertial particles. The 
neutrally buoyant particles cover a Stokes number 
from 0.38 to 1.23, while the heavy particles cover 
a Stokes number from 0.25 to 1.04 (see Table |l|. 
These two ranges overlap allowing some compar- 
ison. As for each class of particles the diameter 
and density are kept constant, the Stokes number 
is varied by tuning the fiow dissipation time scale. 
Therefore, it cannot be varied independently of the 
Reynolds number of the carrier fiow. 

We identify the particles on the images as local 
maxima with light intensity higher than a thresh- 
old, assuming in a first approximation that all the 
particles illuminated in the laser sheet belong to 
one plane. The center of the particles is deter- 
mined as the center of mass of all the pixels sur- 
rounding one local maximum. Due to the high 
contrast between the light diffused by the particles 
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Figure 2: (a) Typical raw acquired image of neutrally buoyant particles and (b) the detected particles located 
in real space with the associated Voronoi diagram (c) Heavy particles located in real space with the associated 

Voronoi diagram (corresponding raw image not shown). 



and the background, changing slightly the thresh- 
old value does not significantly impact the number 
of detected particles, which is of order 100 (resp. 
150) for the neutrally buoyant (resp. heavy) parti- 
cles (note that dealing with finite-size particles, the 
maximum authorized seeding density is drastically 
reduced by particles in the bulk blocking/eclipsing 
the image of particles in the laser sheet, com- 
pared e.g. to experiments with small particles 
where thousands of particles per image are typi- 
cally recorded [15 ). The number of detected neu- 
trally buoyant particles remains constant in time, 
indicating a good stationarity of seeding concentra- 
tion as expected for non-settling particles. How- 
ever large heavy particles tend to settle for low 
impeller rotation rates /imp- The heavy particles 
we consider here are sufficiently large to be con- 
sidered as finite size and sufficiently small to pre- 
vent significant settling as the entrainment by the 
fiow is still capable to keep them in suspension. 
Because of this limitation we did not consider big- 
ger particles, and we did not investigate regimes 
where /imp < 2 Hz (for which settling becomes im- 
portant). Moreover, we make sure that the flow 
is already set in motion when the particles are in- 
serted in the vessel to prevent them from settling 
immediately. Additionally, the number of parti- 
cles per image is monitored, and experiments are 
repeated (after reloading particles) if too many are 
found to have settled (a typical experiment can run 
a few hours with relatively stationary seeding con- 
ditions) . 

The particle concentration field is investigated 
using Voronoi diagrams; this technique recently in- 
troduced for the investigation of preferential con- 
centration of small water droplets in a turbulent 
airfiow |15| was shown to be particularly efficient 
and robust to diagnose and quantify clustering 
phenomena. A given raw image, the detected par- 
ticles and the associated Voronoi diagram are pro- 
vided for neutrally buoyant particles in Fig. |2ja,b) ; 



Fig. [2f^c) shows a typical Voronoi diagram for the 
heavy particle case. The Voronoi diagrams give a 
tesselation of a two-dimensional space where each 
cell of the tesselation is linked to a detected par- 
ticle, with all points of one cell closer to its asso- 
ciated particle than to any other particle. Thus, 
the area of each Voronoi cell is the inverse of the 
local concentration of particles, i.e. Voronoi area 
fields are a measure of the local concentration fields 
at interparticle length scale. To compare the re- 
sults of experiments made with different amount 
of detected particles per image, the Voronoi area 
is normalized using the average Voronoi area A de- 
fined as the mean particles concentration inverse, 
independent of the spatial organization of the par- 
ticles. Therefore, we focus in the rest of the let- 
ter on the distribution of the normalized Voronoi 
area V = A/ A. Clustering properties are quanti- 
fied by comparing the probability density function 
(PDF) of Voronoi cell areas obtained from the ex- 
periments to that of a synthetic random Poisson 
process (RPP) whose shape is well approximated 
by a Gamma distribution |18j . 

The PDFs of Voronoi cell areas for the differ- 
ent experiments described in Tableware plotted in 
Fig. [3] as well as the Gamma distribution approxi- 
mation for a RPP. In the case of neutrally buoyant 
particles, all PDFs collapse and no Stokes number 
dependency is found. An important finding of the 
present work is that these PDFs are almost undis- 
tinguishable from the RPP distribution, mean- 
ing large neutrally buoyant particles do not ex- 
hibit any preferential concentration whatever their 
Stokes number. In the case of heavy particles, 
the PDFs clearly depart from the RPP distribu- 
tion, with higher probability of finding depleted 
regions (large Voronoi areas) and concentrated re- 
gions (small Voronoi areas), which is the signa- 
ture of preferential concentration. Furthermore, 
the shape of the PDF clearly depends on experi- 
mental parameters {St\^ and/or R\). Interestingly, 
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this dependence is stronger for the small Voronoi 
areas tails, whereas the tail for large Voronoi areas 
(corresponding to depleted regions) appears to be 
more robust. This was also observed for small in- 
ertial particles [T^ . The level of clustering can 
be further quantified using the standard deviation 
of the normalized Voronoi areas a\> = \J (V^) — 1, 
plotted in Fig. |4j For neutrally buoyant particles 
we find a constant value cry — 0.53, which is the ex- 
pected value for a RPP. For heavy particles we find 
(Tv > 0.53 for all the experimental configurations 
investigated, revealing the presence of clustering. 
This result is in agreement with previous measure- 
ments |5j that find clustering for large ((/> ~ 4) 
and heavy (F = 1.4) particles. We find the clus- 
tering level to globally decrease as St\^ and/or R\ 
increase, with no hint of maximum clustering for 
around unity. This is contrary to common ob- 
servations: for numerical Stokesian models |13|[H] 
a maximum is found for St ~ 0.6 ; for experi- 
ments with small particles [T51 [1^1 a peak is ob- 
served for iS'i > 1 while [20' found a mild maxi- 
mum for S't ~ 1. If a maximum of clustering ex- 
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Figure 3: (color online) Superposition of the 
normalized Voronoi' area PDFs for six experiments 

with varying Stokes number (plain lines) and a 
Gamma distribution (dashed line). Inserts represent 

the same results with linear ordinates. 
(a) Neutrally buoyant particles (b) Heavy particles. 
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Figure 4: Standard deviations of the normalized 
Voronoi areas versus the Stokes number with error 

bars (plain lines) for V: neutrally buoyant and o: 
heavy particles, to compare with cry^^ ~ 0.53 in the 
case of a RPP. 



ists in our case (which is reasonable assuming that 
tracer behavior is to be recovered for S'th — > 0), 
the peak would be at S'th < 0.25. However, the 
limit S'th {i.e. i?^ — >■ 0) could not be ex- 
plored here due to the settling effects at low R\. 
The clustering properties (Stokes number depen- 
dence and clusters geometry) for such finite-size 
and heavy particles go beyond the scope of the 
present research and will be investigated in future 
experiments. 

Two important conclusions can be drawn from 
these results, (i) While inertial Stokesian mod- 
els predict clustering within the explored range 
of St, this is not observed in the specific case 
of finite-size neutrally buoyant particles. Conse- 
quently, despite they do not behave as tracers, 
such particles are clearly not of the inertial class. 
The absence of clustering also supports experimen- 
tal results on the dynamics of finite-size neutrally 
buoyant particles |71 |SJ [T7| suggesting that simple 
time-response effects are not sufficient to describe 
the particle/flow interaction and that other mech- 
anisms (such as the role of pressure increments 
at the scale of the particle) are to be accounted, 
(ii) Subsequently, the Stokes number by itself can- 
not be taken as sufficient to characterize cluster- 
ing as we have shown that particles with similar 
Stokes number may or not exhibit preferential con- 
centration. Note that the limitation of the Stokes 
number has previously been shown for dynamical 
properties (rather than spatial distribution prop- 
erties as illustrated here) of both isolated parti- 
cles in turbulent flows j31 HI] and two-way cou- 
pling effects [22J. These observations combined to 
ours stimulate the need for further investigations 
on the possible connection between dynamical fea- 
tures and preferential sampling of particles (in- 
cluding for instance turbophoresis and stick-sweep 
mechanisms , but also ergodic mechanisms [23] ) , 
by coupling Voronoi analysis of particles distribu- 



5 



tion to Lagrangian tracking [U [TS] . Acknowledgments 

We currently investigate which set of parameters 
f{(j),T,R\,...) drives the finite-size particle clus- 
tering, in particular by exploring how clustering 

is affected when Stokes number is varied (using This work is partially supported by the ANR 
particles with different sizes) at a fixed Reynolds project TEC. We thank D. Le Tourneau for the 
number. manufacturing of the LEM. 



R. Monchaux, M. Bourgoin, and A. Cartellier, [13 
Int. J. Multiphase Flow 40, 1 (2012). 
S. Goto and J. C. Vassilicos, Phys. Rev. Lett. 100, 
054503 (2008). [14 

E. Calzavarini, M. Kerscher, D. Lohse, and 

F. Toschi, J. Fluid Mech. 607, 13 (2008). [15 
Y. Tagawa, J. M. Mercado, V. N. Prakash, 

E. Calzavarini, C. Sun, and D. Lohse, J. Fluid [16 
Mech. 693, 201 (2012). 

M. Guala, A. Liberzon, K. Hoyer, W. Kinzelbach, 

and A. Tsinober, J. Turb. 9, N34 (2008). [17] 

G. A. Voth, A. La Porta, A. M. Crawford, 

J. Alexander, and E. Bodenschatz, J. Fluid Mech. [18 
469, 121 (2002). 

N. M. Qureshi, M. Bourgoin, C. Baudet, [19 
A. Cartellier, and Y. Gagne, Phys. Rev. Lett. 99, 
184502 (2007). 

R. D. Brown, Z. Warhaft, and G. A. Voth, [20 
Phys. Rev. Lett. 103, 194501 (2009). 

E. Calzavarini, R. Volk, M. Bourgoin, E. Leveque, [21 
J.-F. Pinton, and F. Toschi, J. Fluid Mech. 630, 

179 (2009). 

H. Homann and J. Bee, J. Fluid Mech. 651, 81 [22 
(2010). 

F. G. Schmitt and L. Seuront, J. Marine Systems [23 
70, 263 (2008). 

H. Xu and E. Bodenschatz, Phys. D 237, 2095 
(2008). 



J. Bee, L. Biferale, M. Cencini, A. Lanotte, 
S. Musacchio, and F. Toschi, Phys. Rev. Lett. 98, 
084502 (2007). 

S. W. Coleman and J. C. Vassilicos, Phys. Fluids 
21, 113301 (2009). 

R. Monchaux, M. Bourgoin, and A. Cartellier, 
Phys. Fluids 22, 103304 (2010). 
R. Zimmermann, H. Xu, Y. Gasteuil, M. Bour- 
goin, R. Volk, J.-F. Pinton, and E. Bodenschatz, 
Rev. Sci. Instru. 81, 055112 (2010). 
R. Volk, E. Calzavarini, E. Leveque, and J.-F. Pin- 
ton, J. Fluid Mech. 668, 223 (2011). 
J.-S. Ferenc and Z. Neda, Physica A 385, 518 
(2007). 

M. Obligado, M. Missaoui, R. Monchaux, 
A. Cartellier, and M. Bourgoin, J. Phys.: 
Conf. Ser. 318, 052015 (2011). 
J. R. Fessler, J. D. Kulick, and J. K. Eaton, 
Phys. Fluids 6, 3742 (1994). 

N. M. Qureshi, U. Arrieta, C. Baudet, A. Cartel- 
her, Y. Gagne, and M. Bourgoin, Eur. Phys. J. B 
66, 531 (2008). 

F. Lucci, A. Ferrante, and S. Elghobashi, 
Phys. Fluids 23, 025101 (2011). 
K. Gustavsson and B. Mehlig, Europhys. Lett. 96, 
60012 (2011). 



